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We perform a systematic study of the masses of charmed and bottom baryons in the framework 
of the QCD sum rule approach. Contributions of the operators up to dimension six are included 
in operator product expansion. The resulting heavy baryon masses from the calculations are well 
consistent with the experimental values, and predictions to the spectroscopy of the unobserved 
bottom baryons are also presented. 
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I. INTRODUCTION 



During the past several years there has been tremendous progress in the experimental investigations of 
the heavy baryon spectroscopy. With the precise measurement for the mass of A{, by the CDF collaboration 
[l| , the DO collaboration proclaimed the observation of 2f, Q , which was quickly confirmed by the CDF 



collaboration [3]. The first observations of and have been reported by the CDF [J]. The BABAR 
collaboration announced the observation of fi* [sf and the production of fie from B decays M . The Sc as 
well as the excited states of Sc were set forth by the Belle and the BABAR collaborations [71- Therefore, 
a large amount of experimental data on charmed and bottom baryons has become available. 

On the other hand, various theoretical models have been used to study heavy baryon masses, such as 
quark models (sl, Isl, mass formulas [10], and lattice QCD calculations 11]. From QCD sum rules 12). 
masses of the heavy baryons were primarily calculated in the heavy-quark limit fisl , and subsequently in 



the heavy quark effective theory IJ, [la, ll6[ . In Refs. [l7|, the calculations for the heavy baryons began 



with the full theory. Recently the masses of and fig were tested in QCD sum rules 
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renewed works were inspired by the current significant observations continually 
from the motivation to evaluate the spectroscopy of the heavy baryons systematically in full QCD, we 
shall study mass sum rules for the heavy baryons, with the technique developed in 23|, |2J]. For the cases 



of Sq and flr> , the interpolating currents utilized in this work are not entirely the same as the ones used 
in Refs 
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The paper is organized as follows. In SecHIl QCD sum rules for the heavy baryons are derived. Section 
ml contains numerical analysis, a brief summary, and some discussions. 

II. QCD SUM RULES FOR THE HEAVY BARYONS 



It is of interest to apply the QCD sum rule to study the heavy baryons composed of a heavy quark 
{Q — c, b) and two light (u, d, or s) quarks. The basic point is to choose the suitable interpolating current. 
For the ground states, the currents are correlated with the spin-parity quantum numbers J^' = 0"*" and 
J^'^ — 1+ for the light diquark system, along with the heavy quark forming the state with — ^ and 
the pair of degenerate states with — 5^ ^^'^ ^ i^' which may determine the choice of F^ and F^ 
matrices in baryonic currents [3]- For the F^ matrices for the excited baryons with I{J^) — 0(i ), it 
might be referred to the current of the heavy-light vector meson fi^] . For the baryons with J = | , the 
currents may be gained from those of baryons with J = ^ using SU(3) symmetry relations 26[. Thus, the 
following forms of currents are adopted 



17, 25 



jAiQ = £abc{qiaCTkq2b)Tf,Qc 
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JA-Q = Sabc[—^{qiaCTkQb)Tkq2c+--^iqiaCTkq2b)Ti,Qc], 

jsQ = £abc{qL(^^kSb)'r'kQc, 
£abc{qiaC^^kSb)^kQc, 



JeL' = eabc[^{qLcrkQb)r'kSc + ^{qLcrkSb)r'M (i) 



jHiQ = eafcc(giaCTfcSfc)rj.Qc, 

JH^^ = Sabc[^{qLCTkQb)r',s, + ^iqlCTkSb)T'M 

jUq = £abciSaCTkSb)TkQc, 

2 '1 ' 

jn^ = eabc[-^{slCTkQb)TkSc + -j=isJCTkSb)rkQc]. 



Here the index T means matrix transposition, C is the charge conjugation matrix, a, 6, c are color indices, 
and gi, g2 are u or d depending on the concrete quark contents of the corresponding heavy baryons. The 
choice of and matrices are shown in TABLE HI 



TABLE L The choice of and F^ matrices in baryonic currents 
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The QCD sum rules for the heavy baryons are constructed from the two-point correlation function 

n(q2) ^ J I dW'^-(0|T[j-(^)J(0)]|0). (2) 
Lorentz covariance implies that the two-point correlation function has the form 

U{q')^n,{q')+4U2{q^). (3) 

For each invariant function FEi and II2, a sum rule can be obtained, which is shown below. Phenomeno- 
logically, the correlator can be expressed as a dispersion integral over a physical spectral function 

n q2 ^ x^l^ £L + _ / ds ^ + subtractions, (4 

Mfj -q^ IT s - 



where Mh denotes the mass of the heavy baryon. In obtaining the above expression, the Dirac and 
Rarita-Schwinger spinor sums 

Y,N{q.s)N(q,s) = i + MH, (5) 

S 

for spin-^ baryon, and 

N,{q, s)NM, s) = {i + Mh){9,. - I7.7. + ""'^llil^^' - (6) 



H 



for spin- 1 baryon have been used. In the operator product expansion (OPE) side, short-distance effects are 
given by Wilson coefHcients, while long-distance confinement effects are attributed to power corrections 
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and parameterized in terms of vacuum condensates. Hence 

n,(g2) = Il^\q^)+Ilf^^%%^^l,2. (7) 

We work at leading order in ag and consider condensates up to dimension six. The strange quark is dealt 
as a light one and the diagrams are considered up to order m^. To keep the heavy-quark mass finite, 
the momentum-space expression for the heavy-quark propagator is used. We follow Refs. [2^, 24 1 and 



calculate the light-quark part of the correlation function in the coordinate space, which is then Fourier- 
transformed to the momentum space in D dimension. The resulting light-quark part is combined with the 
heavy-quark part before it is dimensionally regularized at D = 4. For the heavy-quark propagator with 



two and three gluons attached, the momentum-space expressions given in Ref. [25j are used. With Eq, 



([7]), the correlation function in the OPE side in terms of a dispersion relation can be written as 

J ml s-q^ 

where the spectral density is given by the imaginary part of the correlation function 

p,{s) = ilmnpPE(,). (9) 

After equating the two expressions for n(q^), assuming quark-hadron duality, and making a Borel trans- 
form, the sum rules can be written as 

\]jMHe-^'»/^'" = r dspi{s)e-'/^'" +Bn'i''^'^, (10) 

To eliminate the baryon coupling constant Xh and extract the Mh, first take the derivative of Eq. pU]) 
with respect to divide the result by Eq. pI7|) itself, and similarly deal with Eq. pTjl to yield 

Mf, = {^ysp,is)se-/'''+d/di^^)Mr''Hs)}/{^ysp^is)e-^/'^^ (12) 
M], = { [^\sp2is)se-^/'''' +d/di-^)M'i'''^'^{s)}/{ [^\sp2is)e-^^^^^^^ (13) 

where 



with 



■Q "'Q 



pP^^*(s) = ^l^ruQ J^da{^—^ f{ml- saf, (15) 
pf\s) = ^mj\a[i'—^r + 2], (16) 

Z TT J (X 
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^ ™ ' da\\(\ + 2a) - ^]e-"Q/M^ )^ (90) 



pP^'^^s) = -_-L_^ da{^—^ f{l - 3a){mQ - saf - ^msniQ da ^^ (mg - sa), (91) 

pr^(^) - 3^™. rda(3-5a)-^mQ(l-!^)2, (92) 
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III. NUMERICAL RESULTS AND DISCUSSIONS 

In the numerical analysis, the input values are taken as = 1.25±0.09 GeV, to^ = 4.20±0.07 GeV [2^ 
with rus = 0.13 GeV, {qq) = -(0.23)3 GeV^ (ss) = Q.^qq), (gqcr-Gq) = TOg(gg),TOg = 0.8 GeV^ (52G2) = 



0.5 GeV^ and (.9303) 0.045 GeV^ [24|. The proper ranges of the thresholds can be determined, to 
which the stability of the Borel curves should not be sensitive. The Borel windows are fixed in this 
way: the lower limit constraint for A/2 is obtained from the condition that the perturbative contribution 
should be larger than the condensate contributions while the upper one is got by the pole contribution 
larger than the contimmm contribution 29]. Giving an illustration, the comparison between pole and 
continuum contributions from Eq. (12) for ^/sq = 3.4 GeV for ilc is shown in Fig. 1(a), and its OPE 
convergence by comparing the different contributions is shown in Fig. 1(b). The analysis for the others 
has similarly been done, but the corresponding figures are not listed to keep the paper from being too 
lengthy. Accordingly, the threshold are taken as those values exhibited in corresponding figures, and Borel 
windows are M'^ = 1.5 ~ 3.0 GeV^ and 4.5 ^ 6.0 GeV^ for charmed and bottom baryons, respectively. 
The Borel curves for the dependence on Af2 of the heavy baryon masses are shown in Figs. 2-6. In Table 
im we present our results for the masses of the heavy baryons and compare with experimental data and 
other theoretical approaches. In the numerical evaluation, the results gained from the two sum rules have 
been averaged to decrease the systematic errors. It is worth noting that the uncertainty in the results are 
merely owing to Borel windows, not involving the ones rooting in the variation of the quark masses and 
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QCD parameters. 

In conclusion, we have employed the QCD sum rule approach to calculate the masses of charmed and 
bottom baryons including the contributions of the operators up to dimension six in OPE. The final results 
extracted from the sum rules are well compatible with the existing experimental data. Predictions to the 
spectroscopy of the unobserved bottom baryons arc also presented. Although there has been enormous 
progress in experimental aspects and many theoretical works have been done for the heavy baryons, plenty 
of problems are desiderated to resolve. It is worth elucidating that most of the quantum numbers for 
the heavy baryons have not been determined experimentally, but are assigned by PDG on the basis of 
quark model predictions, which arc waiting for experimental identification, especially for several higher 
excited states. More data on bottom baryons are earnestly expected by putting into operation the Large 
Hadron Collider, which may supply a gap of experimental data in the near future. From the theoretical 
point of view, it might be meaningful to reanalyze the QCD sum rules in full theory for the heavy baryons, 
taking into account the QCD 0(0;^) corrections to improve the results, which are not involved in this work. 
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FIG. 1: In (a), the dashed line shows the relative pole eontribution (the pole contribution divided by the total, 
pole plus continuum contribution) and the solid line shows the relative continuum contribution from Eq. (12) for 
y/so = 3.4 GeV for Q,c. Its OPE convergence is shown in (b) by comparing the perturbative, quark condensate, 
two-gluon condensate, three-gluon condensate, and mixed condensate contributions. 
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FIG. 2: The dependence on M for the masses of Aic, A^^, Ai;,, and A^^. The continuum thresholds are orderly 
taken as ^ = 3.1 ~ 3.3 GeV, ^ = 3.2 - 3.4 GeV, ^ = 6.5 - 6.7 GeV, and ^ = 6.6 - 6.8 GeV. (a) and 
(c) are from the sum rule (|12p . (b) and (d) from the sum rule (|13|l . 





FIG. 4: The dependence on for the masses of S^, ^ ^fe^ a-nd S^,*. The continuum thresholds are orderly taken 
as ^ = 3.1 - 3.3 GeV, ^ = 3.2 - 3.4 GeV, ^ = 6.6 - 6.8 GeV, and ^ = 6.6 - 6.8 GeV. (a) and (c) are 
from the sum rule (|12|l . (b) and (d) from the sum rule (|13p . 
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FIG. 5: The dependence on for the masses of Sic, ^ic^ —ib? ^ih- The continuum thresholds are orderly 
taken as ^ = 3.4 ~ 3.6 GeV, ^ = 3.5 ~ 3.7 GeV, ^ = 6.7 - 6.9 GeV, and ^ = 6.8 - 7.0 GeV. (a) and 
(c) are from the sum rule (|12p . (b) and (d) from the sum rule (|13|l . 
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FIG. 6: The dependence on for the masses of ilc, f^ci ^foi ^h- The continuum thresholds are orderly taken 
as ^ = 3.3 - 3.5 GeV, = 3.4 - 3.6 GeV, ^ = 6.6 - 6.8 GeV, and ^ = 6.7 - 6.9 GeV. (a) and (c) are 
from the sum rule (112^ . (b) and (d) from the sum rule (|13p . 
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TABLE II: The mass spectra of charmed and bottom baryons (mass in unit of MeV except for "Our works") 
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